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Modelling correlations in heavy ion collisions is reviewed from model independent 
characterization of the two-particle correlations to Buda-Lund hydro model fits to 
correlations and spectra in Au + Au collisions at y/s = 130 AGeV at RHIC. 


1 Recent results in modelling correlations 

Modelling correlations in heavy ion reactions is a broad and rapidly expanding 
field that is reviewed here very briefly, due to space limitations. For a more 
detailed and illustrated version of this talk, see ref.El. For more detailed recent 
reviews see also refsJ3’0. 

Although some claims were made in the literature that the shape of the 
two-particle Bose-Einstein co&Kelation functions has to be a (multi-variate) 
Gaussian, by now it is clearoB that these claims were mistaken and even 
approximately Gaussian correlation functions may contain subleading non- 
Gaussian corrections, for example oscillators prefactorsB that are particularly 
important in case of effective binary sourcesH. It is also clear that refined anal¬ 
ysis of the two-particle correlation functions indicates significaid deviations 
from Gaussian shapes in elementary particle induced reactionaj. Sensitive 
searches for non-Gaussian components have not yet been performed in cor¬ 
relation measurements in heavy ion collisions at GERN SPS and at RHIG, 
to our best knowledge. A model independent parameterization of the two- 
particle copcfilation function is given in terms ^ the Edgeworth or Laguerre 
expansionsEB. These expansions were shownH’B’Ed to be sensitive to even small 
deviations from Gaussian or exponential shapes. 

The deviation from an exactly Gaussian form of the two-particle correla¬ 
tion function can be caused by the decay products of long-lived resonanc^. 
A core-halo picture has been introduced to describe such a situationD’LllS’Q. 


csorgo-datong: submitted to World Scientific on February 1, 2008 1 



Within the core-halo picture, the intercept parameter A* of the two-particle 
correlation function is interpreted as the squared, momentum dependent frac¬ 
tion of particles emitted from the core of the interaction, and the shape of 
the correlation function is interpreted as a measure of the core part of the 
interaction region. The core-halo structure of the particle emitting source has 
recently been observed with the help of a new imaging method, developed by 
D. A. Brown and P. Danielewitaij, that reconstructs the relative source dis¬ 
tribution directly from the measured data points with the help of the known 
final state interactions of the emitted pair. 

Good progress has been made recently in the study of multi-particle 
Coulomb corrections. A wave-function integration method has been devel¬ 
oped for S^iparticles, utilizing a cluster decomposition of the multiparticle 
final statell3ll3. Such Couloii^corrections have to be applied within the con¬ 
text of a core-halo picturaiil’Ej. A modification of the Coulomb interactions 
(permutations in the charge allocation) has recently been proposed to develop 
Monte-Carlo event generators with Bose-Einstein correlationsEEl. 


An interesting new direction of modelling correlations in hea^zy. 


sions is the study of b. 
well as fermionic cases 


;-to-back correlations (BBC) for bosonicO’l 


dli- 
j as 


(bBBC and fBBC). When the thermalized vacuum 
state decays, it always produces particle - anti-particle pairs with opposite 
spin and momenta. Hence the strength of BBC is increasing inversely j-ssfith 
the single particle spectra and in principle, it can be unlimitedly largal3’EJ. 
It would be certainly of great interest to observe BBC experimentally. The 
bosonic BBC in heavy ion physics has a famous aflakigy in astropysics, namely 
the Bekenstein-Hawking radiation of black holesc3’E3. Both effects are related 
to squeezing and decay of a modified vacuum to pairs of bosons of asymptotic 
fields with back-to-back momentum and opposite quantum numbers. 


2 Hydro models fit spectra at RHIC 

Modelling high energy heavy ion collisions is not only difficult because of the 
large number of degrees of freedom involved in the process but also because 
various theoretical concepts and approaches are relevant for different stages of 
the collision. _A comprehensive review of this process has recently been given 
by S.A. Bas^J. During a partonic cascading process (local) thermalization 
is achieved, that can be maintained for a while due to the intensive collision 
rate and the time evolution corresponds to nuclear fluid dynamics. As the 
system rarities, the collision rate within the “macroscopically infinitesimal” 
fluid cells decreases and one expects that the hydrodynamical approximation 
may break down and non-equilibrium hadronic transport models start to play 


csorgo-datong: submitted to World Scientific on February 1, 2008 2 





aroleB. A version of this well established concept has been utilized to describe 
the single particle spectra and two-particle correlations in 130 AGeV An -I- 
Au collisions at RHICcJ. Numerical solution of relativistic hydrodynamics 
was terminated by particle freeze-out and a subsequent resonance decay and 
hadronic cascading by a code developed by Soff, Bass and Dumitru. The 
model fitted the singlcjuarticle spectra measured by the PHENIX and STAR 
experiments at RHIC^a, but it overestimated, in a statistically unacceptable 
manner, the experimentally observed lon^udinal and out radius components 
of the Bose-Einstein correlation functionsaa. Both the longitudinal and the 
out radius components are sensitive to the distribution of particle production 


. with the calculation could be its utilization of 
of the source distribution_tha±_are known to 


in time. One of the prc 
the Gaussian variances_ 

overestimate the HBT radii for core-halo type of systemsBuc 

An interesting study of the identified jparticle elliptic flow at RHIC has 
been reported by the STAR collaboratiorO. Erom the mass dependence of 
the second flow coefficient as a function of transverse momentum, a sugges¬ 
tive evidence has been obtained that early local equilibration followed by a 
hydrodynamic expansion is responsible for the particle spectra in central and 
mid-peripheral collisions up to pt < 1.5 GeV. A similar conclusion has been 
reached by Kolb, Huovinen, Heinz and Heiselberg, who observed in addition 
that the very rapid thermalization in Au -I- Au collisions at RHIC provides 
a serious chal^ge for kinetic approaches based on classical scattering of on- 
shell particlesES. On the other hand, the high pt behaviour of the elliptic 
flow pattern was shown to be a sensitive measure of the initial parton den¬ 
sity distribution of an initial quark-gluon plasma phaseEil. Detailed numerical 
solution of relativistic hydrodynamics has been reported to describe well the 
single particle spectra, the first and the second J.ow coefficients for various 
domains of centrality and transverse momentumE3, however, this calculation 
also over-predicted by a factor of 2 the measured, effective HBT radii of the 
two-particle Bose-Einstein correlation function, possibly indicating a problem 
with the freeze-out time distribution in this kind of calculations. A good de¬ 
scription of various aspects of the measured single particle spectra, and the 
rrit dependence of the longitudinal HBT xadius parameters has been reported 
by Hirano, Morita, Muroya and Nonaka^ based on an exact solution of rel¬ 
ativistic hydrodynamics. The order of magnitute of the calculated side and 
out radii were also found to be in agreement with the experimental data, 
however, the side radius component was slightly underestimated and the out 
radius component over-estimated by this model, and the observed decrease of 
the out component with increasing transverse mass has not been reproduced. 
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Figure 1. Simultaneous fit of the Buda-Lund hydro model to preliminary data of 
PHENIX^®’^^ (left panels), and (right panels) on particle spectra and corre¬ 

lations of central Au -|- Au collisions at RHIC, \/s = 130AGeV. 


3 Buda-Lund hydro fits spectra and HBT radii at RHIC 

It has been observed that simple parameterizations of the freeze-out phase 
space distribution of a locally thermalized, three dimensionally expanding, 
cylindrically symmetric finite sy^ms describe the preliminary single particle 
spectra at RHIC reasonably wella in the whole available kinematic domain. 

The Buda-Lipd hydrodynamical parameterization (BL-H) has been de¬ 
veloped in refs.E§E3 to describe particle correlations and spectra in heavy ion 
collisions at CERN SPS energies. The model was formulated on the principle 
of cylindrical symmetry in central collisions and a relativistic, 3 dimensional 
flow profile that corresponds to a Hubble type of transverse flow sitting on top 
of a longitudinal Bjorken flow. Under certain conditionsES, the BL-H param¬ 
eterization predicted an i?; ~ ~ i?o oc l/ydrij scaling of the HBT radius 

parameters, as a direct reflection of the cylindrical symmetry of the source. 

Both the single particle spectra and thp tw o -narticle Bose-Einstein corre¬ 
lation functions are analytically calculated l3’t§EZliE3 as a function of 8 model 
parameters plus constants of normalization for the absolutely normalized sin¬ 
gle particle spectra. The BL-H model fits simultaneously the double differen¬ 
tial single particle spectra and the mt dependence of the radius parameters 
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Table 1. Source parameters from simultaneous fittings of CERN SPS (NA49, NA44 and 
preliminary WA98) and RHIC (preliminary PHENIX and STAR) particle spectra and HBT 
radius parameters with the Buda-Lund hydrodynamical model. 


BL-H 

parameters 

STAR prel. 

PHENIX prel. 

{RHIC) 

{SPS) 

Value 

Error 

Value 

Error 

Value 

Error 

Value 

Error 

To [MeV] 

143 

± 

4 

140 

± 

3 

142 

± 

2 

139 

± 

6 

{ut) 

0.76 

± 

0.06 

0.67 

± 

0.3 

0.71 

± 

0.05 

0.55 

± 

0.06 

Rg [fm] 

7.5 

± 

0.3 

6.6 

± 

0.3 

7.1 

± 

0.7 

7.1 

± 

0.2 

TO [fm/c] 

8.8 

± 

0.5 

7.8 

± 

0.3 

8.3 

± 

0.7 

5.9 

± 

0.6 

At [fm/c] 

0.01 

± 

1.2 

1.3 

± 

1.0 

1.3 

± 

1.0 

1.6 

± 

1.5 

Ar) 

1.0 

± 

0.9 

1.5 

± 

0.1 

1.3 

± 

0.4 

2.1 

± 

0.4 

{^)r 

0.09 

± 

0.02 

0.01 

± 

0.02 

0.05 

± 

0.05 

0.06 

± 

0.05 

{^)t 

1.4 

± 

0.4 

0.9 

± 

0.1 

1.2 

± 

0.4 

0.59 

± 

0.38 

X-‘INDF 

46/58 

= 0.79 

45/54 

= 0.83 

0.81 

1.20 


in hadron-proton reactions at CERN SPSB. The BL-H also describes the 
NA49, NA44 and the preliminary WA98 data for spectr|a_and correlations at 
CERN SPS, with a statistically acceptable fittD, as summarized in 

the last column of Table 1 for a comparision with BL-H fits to the prelimi¬ 
nary PHENIX and STAR data for Au -I- Au at ^/s = 130 AGeV at RHIC. 
We observe that the BL-H fits of these RHIC data are satisfactory. The cen¬ 
tral value of the freeze-out temperature, Tq is about the same at RHIC as at 
CERN SPS, and within three standard deviations, most of the fit parameter 
values are similar. However, the mean value of the freeze-out time tq, and 
the strength of the transverse flow at the geometrical radius, (ut), is signifi¬ 
cantly larger at RHIC than at CERN SPS. The system at RHIC freezes out 
later, with a larger transverse flow parameter (itt) than at CERN SPS. The 
preliminary results indicate that within the presently big errors the trans¬ 
verse geometrical size Rg is about as big at RHIC as at CERN SPS. This 
research was supported by grants OTKA T026435 T034296 N25487 (with 
NWO), NSF-MTA-OTKA 0089462 and US DOE DE - FG02 - 93ER40764. 
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